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ABSTRACT 

CL20,  a  high-energy  explosive  compound,  is  a 
polyazapolycyclic  caged  polynitramine  (2,4,6,8,10,12- 
Hexanitro-  2,4,5,8,10,12-hexaazaisowurtzitane)  [Nielson, 
et  al.  1998],  The  combustion  and  detonation 
characteristics  of  CL20  can  be  improved  if  it  is  formed 
into  nanoparticles  of  uniform  size.  A  new,  promising 
process  for  particulation  of  materials  utilizes 
environmentally  benign  compressed  gases  as  either 
solvents  or  anti-solvents.  Predictive  models  are  required 
to  describe  the  solubility  and  phase  behavior  of 
supercritical  solutions  of  CL20  and  supercritical  carbon 
dioxide  and  for  process  simulation  and  development. 
Here,  the  solubility  of  CL20  in  supercritical  carbon 
dioxide  was  evaluated  using  the  Peng-Robinson  cubic 
equation  of  state.  Critical  properties,  vapor  pressure,  and 
other  required  thermodynamic  properties  were  estimated 
using  a  variety  of  available  estimation  techniques, 
including  the  group  contribution  methods  of  Lydersen  and 
of  Joback.  A  Fortran  program  to  predict  the  solubility  of 
CL20  was  developed  during  the  course  of  this  project. 
The  program  was  validated  using  available  literature  data 
for  the  solubility  of  naphthalene  and  of  biphenyl  in 
supercritical  carbon  dioxide.  The  applicability  of  the 
estimation  techniques  employed  for  the  critical  properties 
for  CL20  was  established  using  these  same  techniques  to 
estimate  the  critical  properties  of  comparable  compounds, 
including  RDX  and  HMX.  Solubility  data  for  RDX  in 
supercritical  carbon  dioxide  reported  in  the  literature  were 
also  used  to  establish  the  validity  of  the  estimation 
approach.  Solubility  was  predicted  over  the  temperature 
range  of  305.15  to  368.15  K  and  over  the  pressure  range 
of  74  to  150  atm.  In  general,  as  the  temperature  increases, 
the  solubility  decreases,  while  as  the  pressure  increases, 
the  solubility  increases.  For  CL20,  the  estimated  vapor 
pressures  are  extremely  small,  on  the  order  of  10"18  at 


ambient  temperature,  increasing  to  10"13  at  368.15  K. 
Thus,  the  predicted  solubilities  are  also  small  (range  of 
10' 13  to  10"6),  with  the  highest  solubility  predicted  for 
308.15  K  (35°C)  and  150  atm.  In  the  region  of 
temperature  nearest  the  critical  point  of  carbon  dioxide, 
the  influence  of  the  supercritical  fluid  is  stronger.  It 
would  be  most  desirable  to  operate  in  this  region  if  one  is 
interested  in  maximizing  the  solubility  of  the  solute, 
CL20,  in  the  carbon  dioxide. 

1.  INTRODUCTION 

The  Army  is  interested  in  developing  a  green  process 
to  produce  nanocrystalline  (<  300nm)  particles  of  the  high 
energy  explosive  compound  e-CL20  (2,4,6,8,10,12- 
Hexanitro-2,4,5,8,10,12-hexaazaisowurtzitane).  The 
combustion  and  detonation  characteristics  of  the  CL20 
can  be  improved  if  it  is  formed  into  nanoparticles  of 
uniform  size.  Processes  such  as  Supercritical  Antisolvent 
process  (SAS)  and  condensation  with  Rapid  Expansion  of 
Supercritical  Solutions  (RESS),  Gas  Antisolvent  (GAS), 
and  Supercritical  Antisolvent  with  Enhanced  mass 
transfer  (SAS-EM™)  are  potential  green  processes  for 
producing  ultrafine  particles.  In  these  processes,  the 
material  to  be  particulated  will  be  dissolved  (solubilized) 
into  an  environmentally  benign  solvent  such  as 
supercritical  carbon  dioxide  and  then  condensed  to 
ultrafine  particles  by  reducing  the  pressure  and 
temperature  of  the  mixture. 

Theoretical  and/or  predictive  models  are  required  for 
process  simulation  and  to  describe  the  solubility  and 
phase  behavior  of  the  supercritical  solutions  of  CL20  and 
supercritical  carbon  dioxide.  Experimental  data  on  the 
solubility  and  phase  behavior  of  CL20  in  supercritical 
carbon  dioxide  are  not  available.  This  paper  presents  a 
method  to  estimate  the  required  vapor  pressure  data,  other 
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parameters  required  for  the  model,  and  predicted  Table  2.  Thermodynamic  Information  for  HMX  (compiled 
solubility  of  CL20  and  its  close  cousin  RDX.  from  Beilstein  Search). 


2.  LITERATURE  REVIEW 

A  comprehensive  literature  review  was  conducted  to 
identify  sources  of  thermodynamic  information  for  CL20 
and  similar  energetic  compounds,  including  RDX  and 
HMX.  These  energetic  compounds  were  included  since 
they  have  been  used  as  explosive  agents  for  a  longer  time 
and,  thus,  there  would  most  likely  be  more  information 
available.  The  available  data  for  these  compounds  could 
then  be  used  to  validate  estimation  techniques  for  the 
thermodynamic  properties  of  CL20.  Available 
information  on  the  three  compounds  was  compiled  using 
the  CAS  registry  numbers  (CL20,  CAS#  135285-90-4; 
RDX,  CAS#  121-82-4;  HMX,  CAS#  2691-41-0).  Of 
particular  interest  was  the  determination  of  what  data 
were  available  in  the  literature,  including:  the  critical 
point  (temperature,  pressure,  volume,  compressibility), 
the  melting  point,  the  acentric  factor,  the  latent  heat  of 
vaporization,  vapor  pressure,  and  liquid  and/or  solid 
density.  Two  computerized  database  searches  were 
conducted,  yielding  the  following  relevant  information. 


Property 

Value 

Reference 

Critical 

Temperature 

654°C 

Maksimov  (1992) 

Critical  Volume 

611000  cnr'/gmol 

Maksimov  (1992) 

Enthalpy  of 
Vaporization 

158599.97  J/mol  @ 
234°C 

Behrens  (1990) 

Melting  Point 

Range:  198  to 

205°C 

Suri  and  Chapman 
(1988) 

Burov  et  al.  (1999) 

2.2  Chemical  Abstracts  Database 

The  Chemical  Abstracts  Database  revealed  additional 
information  for  the  compounds  of  interest.  For  RDX,  two 
items  of  significance  were  identified.  Dionne  et  al. 
(1986)  measured  the  vapor  pressure  of  RDX  as  a  function 
of  temperature  over  the  temperature  range  37  to  102°C. 
Morris  (1998)  measured  the  solubility  of  RDX  in  dense 
carbon  dioxide  over  the  temperature  range  of  303  to  353 
K.  These  data  provide  a  means  to  assess  the  validity  of 
the  approach  taken  to  estimate  vapor  pressure  and  to 
predict  the  solubility  in  supercritical  carbon  dioxide  for 
the  compound  of  interest. 


2.1  Beilstein  Database 

For  CL20,  the  only  relevant  data  contained  in  the 
database  were  crystal  density  (solid  density)  and  melting 
point.  No  references  were  reported  associated  with  the 
critical  point,  the  normal  boiling  point,  or  phase  change 
properties  (latent  heat,  vapor  pressure).  For  RDX, 
relevant  references  were  found  for  the  critical 
temperature,  critical  volume,  vapor  pressure,  melting 
point  and  enthalpy  of  formation.  The  values  obtained 
from  the  references  are  compiled  in  Table  1.  For  HMX, 
relevant  references  were  also  identified  for  these 
properties,  excepting  vapor  pressure.  The  values  obtained 
from  the  references  are  compiled  in  Table  2. 


Table  1 .  Thermodynamic  Information  for  RDX  (compiled 
from  Beilstein  Search). 


Property 

Value 

Reference 

Critical 

Temperature 

567°C 

Maksimov  (1992) 

Critical 

Volume 

442000  cnr’/gmol 

Maksimov  (1992) 

Vapor 

Pressure 

.000036  -  .0004 

mm  Hg  @  110.6°C 
-138. 5°C 

Edwards  (1953) 

Melting 

Point 

Range:  272°C  to 
278°C 

Burov  et  al. 
(1999) 

3.  THERMODYNAMIC  FRAME  WORK 

The  basis  for  predicting  the  solubility  of  a  solute  in  a 
supercritical  fluid  solvent  is  the  equivalence  of  fugacities 
for  the  particular  solute  in  each  phase: 

/,'=//  (i) 

The  superscript  s  represents  the  solid  phase  and  f  the 
supercritical  fluid  phase.  If  the  solubility  of  the 
supercritical  fluid  in  the  solid  phase  is  assumed 
negligible,  then  the  fugacity  of  the  solute  in  the  solid 

phase,  j\s ,  is  equal  to  the  fugacity  of  the  pure  solute, 

f.  .  The  fugacity  of  the  pure  solute  in  the  solute  phase  is 
evaluated  using  [Sandler  1989]: 

(  fsat  A  T  rP  (  Vs  ^ 

(2) 

The  molar  volume  of  the  pure  solute,  V/ ,  is  assumed 
constant.  The  fugacity  of  the  pure  solute  in  the  saturated 
state,  f-at ,  is  combined  with  the  pressure  in  the 
denominator  to  form  the  fugacity  coefficient  of  the  pure 
solute  in  the  saturated  state,  <f)  'a'  .  Integration  of  equation 
(2)  with  these  substitutions  yields: 
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For  many  low  volatility  compounds,  the  fugacity 
coefficient  for  the  pure  solute  in  the  saturated  state  is  very 
nearly  unity  due  to  the  extremely  low  vapor  pressure  of 
these  compounds  at  ambient  conditions.  The  fugacity  of 
the  solute  in  the  supercritical  fluid  phase  is  given  by: 

//  =  yAp  (4) 


where  y.  is  the  mole  fraction  of  solute  in  the  supercritical 
fluid  phase,  also  defined  as  the  solubility  of  the  solute  in 

the  supercritical  fluid;  (f)j  is  the  fugacity  coefficient  for 


the  solute  in  the  supercritical  fluid  phase;  and  P  is  the 
system  pressure. 


Combination  of  equations  (3)  and  (4)  yields  an  expression 
for  the  solubility  of  the  solute,  y(,  in  the  supercritical  fluid 
[Sandler  1989]: 
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The  ideal  solubility  of  the  solute,  a  function  of 
temperature  and  pressure,  is  represented  by  the  first 
bracketed  term  in  the  expression.  For  the  compounds  of 
interest,  the  ideal  solubility  will  be  extremely  low  due  to 
the  low  vapor  pressure  of  the  solute  and  the  high  pressure 
required  to  achieve  supercritical  conditions.  Non-ideal 
behavior  of  the  supercritical  fluid  phase  is  represented  by 
the  second  bracketed  term  in  the  expression.  The  last 
term  is  the  Poynting  factor,  which  represents  the  effect  of 
pressure  on  the  solid  phase.  The  last  two  bracketed  terms, 
when  combined,  are  known  as  the  enhancement  factor. 
This  factor  represents  the  increase  in  solubility  due  to  the 
solvent’s  supercritical  state. 


The  fugacity  of  the  solute  in  the  supercritical  fluid 
can  be  evaluated  using  a  cubic  equation  of  state  such  as 
the  Redlich-Kwong-Soave  or  the  Peng-Robinson 
equations.  In  this  work,  the  Peng-Robinson  equation  is 
utilized.  The  Peng-Robinson  equation  for  a  mixture  is 
[Peng  and  Robinson  1976;  Modell  and  Reid,  1983]: 


P  = 


RT 


lM 


V  b  M  V(V  +  bM)  +  bM(V-bM) 


(6) 


where  V  is  the  molar  volume  of  the  mixture  while  aVI  and 
bM  are  mixture  dependent  parameters.  For  the  Peng- 
Robinson  equation  of  state,  these  mixture  parameters  are 
evaluated  from  mixing  rules  where  the  pure  component 
analogues,  calculated  from  the  critical  properties,  the 
acentric  factor  ro,  and  the  reduced  temperature  for  each 
species,  are  combined  through  appropriate  combinatorial 
relationships.  The  standard  Peng-Robinson  mixing  rules 
(Equations  6  and  7)  are  utilized  in  this  work. 


N  N 

aM  --  ZZ  V'<> 

;= l  j= l 

hM  =  Z-rA 


(8) 


i= 1 


with: 


a  = 

v 


(l  -  Sy  )^a~aj  and  ati  =  at 


The  binary  interaction  parameter,  8  ,  is  specific  for  the  i- 

u 

j  binary  pair.  In  these  expressions,  the  ai  and  h  represent 
the  pure  component  parameter  values,  which  are 
evaluated  using: 

ai(o>,Tr)  =  a(Tc)a(a>,Tr)  (9) 

0.45724 R2T; 

a<Tc)  = - - - -  (10) 

±  c 

a(a>,Tr) 

=  [l  +  (0.37464  +  1.54226®  -  0.26992® 2  )(l  -  (n) 

=M-VF| 


0.07780 RTc 

Pc 


(12) 


The  Peng-Robinson  equation  is  often  written  in  terms  of 
the  compressibility. 


Z  =  (PV/RT): 

Z3  +Z2(B-l)+z(-3B2 -2B  +  A)+  (13) 
(. B 3  +B 2  -Ab)=0 


where  A  and  B  are  defined  as: 


A=yp„ 

(14) 

(rt) 

buP 

B=  M 

(15) 

RT 


The  use  of  this  equation  of  state  for  evaluation  of  the 
fugacity  coefficient  of  the  solute  in  a  supercritical  fluid  is 
well-documented  in  the  literature  [McHugh  and  Paulaitis 
1980].  Following  the  notation  in  Modell  and  Reid  [1983], 
the  fugacity  coefficient  for  species  i  in  solution  is  found 
through: 


ln(^i ) =  T~~  (Z  - 1)  -  ln(Z  -  B) 


+ 


A 

2V2 B 


In 


Z  +  B  1  -  V2 

- j - < 

Z  +  B  1  +  V2 


(16) 


v  J 

The  prediction  of  solubility  by  this  set  of  equations 
requires  an  iterative  approach.  In  general,  the  temperature 
and  pressure  are  known,  as  well  as  an  estimate  of  the 
solubility,  y;.  If  experimental  data  are  not  available  for 
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use  as  the  initial  estimate,  then  an  estimate  is  calculated 
by  assuming  that  ideal  behavior  is  valid  (i.e., 

yt  =  Plat  IP).  The  estimated  y,  is  then  used  with  the 

given  temperature  to  evaluate  the  pure  component  and 
mixture  parameters.  The  compressibility  is  then 
determined  through  equation  (13)  using  the  given  T  and  P 
with  the  mixture  parameters  through  solution  of  the  cubic 
equation  of  state.  In  this  implementation,  the  cubic 
equation  is  solved  analytically.  The  compressibility  is 
substituted  into  the  expression  for  the  fugacity  coefficient, 
equation  (16),  along  with  the  necessary  pure  component 
and  mixture  parameters.  The  predicted  solubility  is  then 
calculated  through  equation  (5).  The  iterative  procedure 
is  continued  until  the  predicted  solubility  at  the  end  of 
iteration  is  equal  to  the  estimate  at  the  start  of  the  iteration 
to  within  some  prescribed  tolerance. 


4.  ESTIMATION  OF  CRITICAL  AND  OTHER 
REQUIRED  THERMODYNAMIC  PROPERTIES 


Estimation  of  the  critical  properties  and  other 
required  thermodynamic  information  was  undertaken  for 
the  energetic  compounds  RDX,  HMX  and  CL20.  Since  a 
limited  amount  of  pertinent  data  was  available  for  RDX 
and  HMX  (i.e.,  critical  temperature,  critical  volume, 
vapor  pressure),  the  estimation  of  properties  for  these 
compounds  could  serve  as  a  baseline  to  establish  the 
validity  of  the  estimation  procedures  used  and  thus, 
provide  for  greater  confidence  in  the  estimates  obtained 
for  CL20.  Although  CL20  (Figure  1)  is  a  much  more 
complex  molecule  as  compared  to  RDX  (Figure  2)  or 
HMX  (Figure  3),  all  three  contain  -NO2  groups, 
individually  attached  to  a  ring  >N-  group.  Group 
contribution  methods  were  used  to  estimate  the  necessary 
values.  Table  3  shows  the  group  decomposition  for  these 
three  compounds. 
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Fig.  1.  Molecular  Structure  of  CL20 
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Fig.  2.  Molecular  Structure 
of  RDX 


Fig.  3.  Molecular  Structure 
of  HMX 


Table  3.  Group  Decomposition  of  Energetic  Molecules. 


Group  ID 

CL20 

RDX 

HMX 

Number 
of  Type 

Number  of 
Type 

Number  of 
Type 

>N-  (ring) 

6 

3 

4 

-N02 

6 

3 

4 

CH  (ring) 

6 

CH2  (ring) 

3 

4 

The  method  of  Stein  and  Brown  (1994)  was  used  to 
estimate  the  normal  boiling  point  for  each  compound. 
This  method  is  an  extension  of  Joback’s  method  (Joback 
and  Reid  1987)  for  estimation  of  the  normal  boiling  point. 
Stein  and  Brown’s  extension  specifically  involves 
inclusion  of  additional  groups  and  revision  of  Joback  and 
Reid’s  groups  into  smaller  subdivisions.  Of  interest  in 
this  work  is  the  inclusion  of  the  >N-  (ring)  group.  For 
this  group,  the  contributions  to  critical  temperature  and 
pressure  were  identical  to  those  for  a  >N-  (non-ring) 
group,  but  were  different  for  the  normal  boiling  point. 
Stein  and  Brown  use  the  original  relation  of  Joback  and 
Reid,  given  as  equation  (17): 

Tb=  198.2  +  £nlgi  (17) 

i 

An  additional  refinement  to  the  estimate  is  also  given  by 
Stein  and  Brown  and  was  used  in  this  work.  For 
estimated  normal  boiling  points  above  700  K,  the 
corrected  boiling  point  is  given  by: 

Tb  (corr)  =  Tb  -  94.84  +  0.5577Tb  -  0.0007705Tb2  (18) 

Above  700  K,  the  corrected  boiling  point  is  given  by: 

Tb  (corr)  =  Tb  +  282.7  -  0.5209Tb  (19) 

Estimated  normal  boiling  points  calculated  in  this  work 
for  the  compounds  of  interest  are  shown  in  Table  4. 

Table  4.  Estimated  Normal  Boiling  Points 


Compound 

T),  Estimate 
(K) 

Tb(corr) 

Estimate  (K) 

RDX 

717.8 

626.6 

HMX 

891.0 

709.6 

CL20 

1208.7 

861.8 

The  estimation  of  critical  temperature  was  accomplished 
using  Joback’s  method  (Joback  and  Reid  1987)  as 
implemented  in  the  Cranium  software  (Molecular 
Knowledge  Systems,  Inc.).  This  estimation  technique 
requires  the  use  of  the  normal  boiling  point,  and  thus,  the 
values  shown  in  Table  4  were  used.  Estimates  for  the 
critical  temperature  for  the  compounds  of  interest  are 
shown  in  Table  5. 
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The  critical  temperature  estimated  for  RDX,  842.7  K, 
using  the  Tb(corr)  estimate  is  in  excellent  agreement  with 
the  experimental  value  of  567°C  (840  K)  reported  by 
Maksimov  (1992).  Similarly,  the  critical  temperature 
estimated  for  HMX,  913.5  K  using  the  Tb(corr)  estimate, 
is  in  good  agreement  with  the  experimental  value  of 
654°C  (927  K)  reported  by  Maksimov  (1992). 

Table  5.  Estimated  Critical  Temperature  and  Pressure 


Compound 

Critical 

Temperature, Tc 
Estimated  using 
Tb(K) 

Critical 

Temperature, Tc 
Estimated  using 
Tb(corr)  (K) 

Critical 
Pressure, 
Pc  (bar) 

RDX 

965.4 

842.7 

58.0 

HMX 

1147.1 

913.5 

53.0 

CL20 

1483.6 

1057.8 

48.9 

The  estimation  of  critical  pressure  was  also  accomplished 
using  Joback’s  method  (Joback  and  Reid  1987)  as 
implemented  in  the  Cranium  program.  No  literature  data 
for  the  critical  pressure  were  identified  for  the  compounds 
of  interest.  Estimates  for  the  critical  temperature  for  the 
compounds  of  interest  are  also  shown  in  Table  5. 

In  addition  to  the  critical  properties  (temperature  and 
pressure),  the  acentric  factor  is  used  in  describing  the 
temperature  dependence  of  the  attraction  term  in  the 
Peng-Robinson  equation  of  state.  Thus,  estimation  of  the 
acentric  factor  was  also  required  for  the  compounds  of 
interest.  Poling  et  al.  (2000)  recommend  the  estimation  of 
the  acentric  factor  using  the  three  parameter  Pitzer 
expansion: 

ln(Pvp  /Pc)  =  f (0)  +  <of(1)  +  ®2f (2)  (20) 

Neglecting  the  second  order  term,  this  can  be  rearranged, 
as  shown  in  equation  (21)  to  provide  an  estimate  of  the 
acentric  factor  using  the  critical  point  and  the  normal 
boiling  point: 


f'"(Tr) 


f(2)  (Tr) 


- 5.03365t  + 1 . 1 1  505t‘ 5  - 5 .4 12 1  7t2 5  -  7.46628t5  (23) 

T 

- 0.6477 It  +  2.41  539t‘ 5  -4.26979t25  +  3.25259t5  (24) 

T 


In  these  expressions,  t  is  defined  as  (1-Tr).  For  larger 
molecules,  the  inclusion  of  the  second  order  term  can  be 
important.  Thus,  inclusion  of  the  second  order  term  gives 
rise  to  a  quadratic  equation  that  must  be  solved.  This 
quadratic  equation  is  in  the  form: 


orf(2)  +  (of1 (1)  +f(0)  -ln(l. 01325  /Pc)  =  0  (25) 


The  acentric  factor  was  estimated  using  the  values  of 
Tb(corr)  shown  in  Table  4,  and  values  of  Tc  from  Table  5. 
Estimated  values  for  each  of  the  compounds  of  interest 
are  shown  in  Table  6.  These  values  are  compared  to  a 
two  point  extrapolation  of  the  vapor  pressure  curve  to 
obtain  an  estimate  of  the  acentric  factor.  Assuming  a 
linear  relationship  between  logi0(Pr)  and  1/Tr  gives  rise  to 
the  expression  shown  in  equation  (26). 


logio(Pr)  =  -|- 


1-T, 

T 


(26) 


where  B  is  given  by: 


tog  10  (Pc)  -  !°g  lo(P„bp  ) 
1  1 


Using  the  estimated  normal  boiling  point  and  the  critical 
temperature  and  pressure,  the  reduced  vapor  pressure  at  a 
reduced  temperature  of  Tr  =  0.7  was  evaluated  and  then 
used  in  the  definition  of  the  acentric  factor.  Estimated 
values  for  the  acentric  factor  obtained  using  this  method 
are  also  provided  in  Table  6. 

Table  6.  Estimated  Acentric  Factor 


ln(Pc  / 1.01 325) +  f(0,(Tbr) 
fU)(Tbr) 


(21) 


where  Tbr  is  the  reduced  normal  boiling  point  and  pressure 
is  in  bar.  The  functions,  f0)(Tbr),  ^'’(T^),  and  f2)(Tbr)  are 
the  temperature  dependent  correlations  developed  by 
Ambrose  and  Walton  (1989)  for  the  three  parameter 
corresponding  states  method  for  prediction  of  vapor 
pressure.  These  expressions  are  shown  in  equations  (22) 
through  (24)  (Poling  et  al.  2000). 


f«”(Tr) 


-5.97616t  +  1.29874tu  -0.60394t2  5  -1.06841t5 
Tr 


(22) 


Compound 

Acentric 

Factor 

(Estimated 

using 

Cranium) 

Acentric 

Factor 
(Estimated 
using  Eq. 
(22)-(25) 

Acentric 

Factor 
Estimated 
using  Eq. 
(26)-(27) 

RDX 

1.25 

1.22 

1.19 

HMX 

1.57 

1.66 

1.80 

CL20 

2.40 

2.38 

2.18 

The  estimation  of  vapor  pressure  was  also  undertaken  for 
the  model  compounds.  Vapor  pressure  data  were  located 
for  RDX,  but  data  for  HMX  and  for  CL20  were  not 
identified  through  the  literature  search.  Dionne  et  al. 
(1986)  report  the  vapor  pressure  of  RDX  as 
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6473 


(28) 


8.0E-14 


log10(Pvp)  =  22.50- 


where  Pvp  is  in  parts  per  trillion  (v/v)  and  T  is  in  K.  This 
correlation  was  based  on  their  measurements  over  the 
temperature  range  of  37  to  102°C  combined  with 
literature  data  from  three  other  sources.  A  second 
expression  for  the  vapor  pressure  of  RDX  was  (U.S.  Dept, 
of  Labor  2003): 

l°g,.(PJ  =  1414-^  (29) 


where  Pvp  is  in  Pa  and  T  is  in  K.  These  two  expressions 
provide  roughly  the  same  estimates  for  the  vapor  pressure 
of  RDX  over  the  temperature  range  of  interest. 
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Fig.  5.  Estimated  Vapor  Pressure  for  CL20. 


5.  VALIDATION  OF  THE  DEVELOPED  CODE 


Poling  et  al.  (2000)  recommend  the  use  of  the 
corresponding  states  method  with  the  Ambrose-Walton 
expressions  for  the  estimation  of  vapor  pressure.  The 
relevant  expressions  are: 

ln(Pvpr)  =  f(0) +oof<1) +(0f<2)  (30) 

where  Pvpr  is  the  reduced  vapor  pressure,  f is  given  by 
equation  (22),  f n  by  equation  (23),  and  f2)  by  equation 
(24).  Vapor  pressures  were  estimated  according  to 
equation  (30)  using  the  values  for  the  acentric  factor 
given  in  Table  6.  Vapor  pressures  for  RDX  estimated 
using  this  technique  differed  from  the  available  literature 
data  by  a  factor  of  103.  Vapor  pressures  were  also 
estimated  according  to  equation  (26).  The  estimated 
vapor  pressure  is  compared  to  the  available  literature  data 
for  RDX  in  Figure  4.  The  estimates  are  in  good 
agreement  with  the  literature  data  over  the  range  of  25  to 
approximately  60°C.  Above  60°C,  the  estimated  vapor 
pressure  and  the  literature  data  diverge,  with  the  estimated 
vapor  pressure  being  almost  half  that  reported  in  the 
literature  at  80°C.  For  prediction  of  solubility  in 
supercritical  CO2,  it  is  anticipated  that  the  temperature 
range  of  interest  is  between  approximately  35°C  and 
60°C.  Thus,  equation  (26)  was  used  to  estimate  the  vapor 
pressure  of  CL20.  These  estimated  vapor  pressures  were 
then  used  in  the  evaluation  of  the  solubility  of  CL20  in 
supercritical  C02  and  are  plotted  in  Figure  5. 


Fig.  4.  Comparison  of  Estimated  RDX  Vapor  Pressure 
with  Literature  Data. 


The  Fortran  code  was  validated  using  the  solubility 
of  napthalene  and  of  biphenyl  in  supercritical  carbon 
dioxide.  Data  for  these  systems  were  measured  by 
McHugh  and  Paulaitis  (1980).  For  naphthalene  in 
supercritical  C02,  the  value  of  the  interaction  parameter, 
bi2,  was  set  to  0.103,  while  for  biphenyl  in  supercritical 
C02,  a  value  of  0.12  was  used.  These  values  for  the 
interaction  parameter  were  obtained  by  Schmitt  (1984)  by 
regression  of  the  data  of  Tsekhanskaya  et  al.  (1964)  for 
naphthalene,  and  the  data  of  McHugh  and  Paulaitis  for 
biphenyl.  Figure  6  provides  a  comparison  of  the  literature 
data  to  the  predicted  values  for  naphthalene,  while  Figure 
7  provides  the  comparison  of  literature  data  to  predicted 
values  for  biphenyl.  The  predicted  values  are  in  fair 
agreement  with  the  literature  data  for  naphthalene  in  the 
range  of  pressures  from  approximately  the  critical 
pressure  of  carbon  dioxide  (~  74  bar)  to  approximately 
200  bar.  The  predicted  solubilities  at  55°C, 
corresponding  to  the  highest  pressures  examined  by 
McHugh  and  Paulaitis,  diverge  from  the  experimental 
data.  A  reason  for  the  discrepancy  is  that  Schmitt  (1984) 
reported  different  values  of  the  interaction  parameter  for 
each  set  of  isothermal  data.  Thus,  the  same  interaction 
parameter  did  not  provide  the  best  representation  of  data 
for  differing  isotherms. 


Pressure,  P,  (bar) 

Fig.  6.  Solubility  of  Naphthalene  in  Supercritical  C02  - 
Comparison  of  Literature  and  Calculated  Values. 
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Pressure,  P,  (bar) 


Fig.  7.  Solubility  of  Biphenyl  in  Supercritical  C02 
Comparison  of  Literature  Data  and  Calculated 
Values  at  35.8°C 

For  biphenyl  in  supercritical  carbon  dioxide,  the 
agreement  between  predicted  values  is  much  better.  This 
is  most  likely  due  to  the  use  of  an  interaction  parameter 
for  the  data  set  that  had  been  obtained  through  regression 
of  the  same  data  set.  The  primary  motivation  of  using 
these  well-known  data  sets  to  provide  validation  and 
verification  of  the  developed  code  was  that  the 
performance  of  the  code  could  be  more  easily  established 
since  all  required  information  (i.e.,  critical  temperature, 
critical  pressure,  acentric  factor,  vapor  pressure,  saturated 
molar  volume)  was  available  for  the  solutes  (naphthalene 
and  biphenyl)  as  well  as  the  solvent  (carbon  dioxide). 

6.  SOLUBILITY  OF  EXPLOSIVES  IN 
SUPERCRITICAL  CARBON  DIOXIDE 

6.1  Comparison  of  SCF  Predictions  for  RDX  in  C02 
with  Available  Literature  Data 

The  solubility  of  RDX  in  subcritical  and  supercritical 
carbon  dioxide  was  reported  by  Morris  (1998).  The 
reported  values  (in  mg  RDX/g  C02)  were  converted  to 
mole  fractions  for  comparison  with  code  predictions.  The 
initial  estimate  used  in  the  iterative  code  for  the  solubility 
was  the  experimental  value  reported  by  Morris  at  the 
conditions  of  interest.  There  is  little  agreement  between 
the  literature  data  and  the  predictions,  as  shown  in  Figure 
8.  The  literature  solubilities  were  quite  small,  with 
magnitudes  between  10'7  to  10"5,  while  the  predicted 
solubilities  were  larger  with  magnitudes  between  10"7  and 
10"3.  This  is  most  likely  due  to  the  vapor  pressure 
expression  used.  Error  in  the  vapor  pressure  translates 
directly  into  error  in  the  predicted  solubility,  as  evidenced 
through  equation  (5).  If  the  vapor  pressure  is 
overestimated  by  a  factor  of  2,  then  the  predicted 
solubility  will  be  overestimated  to  the  same  degree. 
Accurate  measurement  of  the  vapor  pressure  of  the 


explosive  compound  is  an  essential  requirement  for 
reliable  predictions  of  the  solubility  of  an  explosive  in 
supercritical  carbon  dioxide. 
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Fig.  8.  Solubility  of  RDX  in  Subcritical  and  Supercritical 
Carbon  Dioxide  -  Comparison  of  Predictions  with 
Literature  Data 

6.2  SCF  Predictions  for  CL20  in  Supercritical  C02 

The  solubility  of  CL20  in  supercritical  carbon 
dioxide  was  predicted  using  the  developed  Fortran  code. 
The  estimated  critical  temperature,  critical  pressure  and 
acentric  factor  were  used  in  the  code  along  with  equation 
(26)  for  estimation  of  the  vapor  pressure.  A  two  point  fit 
(critical  point,  normal  boiling  point)  allowed  estimated  of 
the  coefficient  B  in  equation  (26).  Since  no  experimental 
measurements  were  available,  the  initial  estimate  of  the 
solubility  was  set  to  1.0"16.  Solubility  was  predicted  over 
the  temperature  range  of  305.15  to  368.15  K  and  over  the 
pressure  range  of  74  to  150  atm.  The  temperature  range 
corresponds  to  a  range  of  reduced  temperature  of  1.003  to 
1.21  with  respect  to  carbon  dioxide,  while  the  pressure 
range  corresponds  to  a  reduced  pressure  range  of  1.01  to 
2.06.  The  estimated  solubilities  are  shown  in  Figure  9 
(low  temperature)  and  Figure  10  (high  temperature). 
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Fig.  9.  Predicted  Solubility  ofCL20  in  Supercritical  C02  - 
Low  Temperature  Range. 
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Pressure,  P,  (bar) 

Fig.  10.  Predicted  Solubility  of  CL20  in  Supercritical  C02 
-  High  Temperature  Range. 

In  general,  as  the  temperature  increases,  the  solubility 
decreases,  while  as  the  pressure  increases,  the  solubility 
increases.  The  predicted  behavior  is  markedly  different 
near  the  critical  point  of  carbon  dioxide  when  compared 
to  the  behavior  of  naphthalene  in  supercritical  carbon 
dioxide.  For  CL20,  the  estimated  vapor  pressures  are 
extremely  small,  on  the  order  of  10"18  at  ambient 
temperature,  increasing  to  10"13  at  368.15  K.  Thus,  the 
predicted  solubilities  are  also  small  (range  of  10’L'  to  10"6) 
with  the  highest  solubility  predicted  for  308.15  K  (35°C) 
and  150  atm.  In  the  region  of  temperature  nearest  the 
critical  point  of  carbon  dioxide,  the  influence  of  the 
supercritical  fluid  is  stronger.  It  would  be  most  desirable 
to  operate  in  this  region  if  one  is  interested  in  maximizing 
the  solubility  of  the  solute,  CL20,  in  the  carbon  dioxide. 

CONCLUSION 

A  comprehensive  search  of  the  literature  revealed 
limited  data  that  could  be  used  to  validate  the  approach 
taken  in  this  work  for  estimation  of  necessary 
thermodynamic  properties.  The  estimation  of  normal 
boiling  point,  critical  temperature,  critical  pressure, 
acentric  factor,  and  vapor  pressure  for  CL20  was 
accomplished  using  available  group  contribution 
techniques  and  property  definitions.  While  CL20  is  a 
very  complex  species,  comparison  of  estimated  critical 
temperature  for  RDX  and  HMX  to  values  reported  in  the 
literature  provides  a  measure  of  confidence  in  the 
estimated  values  for  CL20.  All  three  species  are 
comprised  of  similar  molecular  groups,  the  primary 
difference  being  the  presence  of  -CH  (ring)  groups  in 
CL20  and  -CH2  (ring)  groups  in  RDX  and  HMX. 

A  Fortran  code  was  developed  for  the  estimation  of 
the  solubility  of  a  solute  in  supercritical  carbon  dioxide. 
Performance  of  the  code  was  validated  using  the  data  of 
McHugh  and  Paulaitis  (1980)  for  naphthalene  and 
biphenyl  solubilities  in  supercritical  carbon  dioxide.  The 
prediction  of  solubility  for  the  explosive,  RDX,  in 
supercritical  carbon  dioxide  was  performed  and  the 


predicted  values  compared  to  available  literature  data. 
The  lack  of  agreement  between  the  predicted  and 
literature  data  is  most  likely  due  to  the  uncertainty 
associated  with  the  vapor  pressure  of  RDX. 

Solubilities  for  CL20  in  supercritical  carbon  dioxide, 
over  a  reduced  temperature  range  of  1.003  to  1.21  and  a 
reduced  pressure  range  of  1.01  to  2.06,  were  predicted 
using  the  developed  code.  These  predictions  lead  to  two 
primary  conclusions.  First,  the  solubility  estimates  are 
very  strongly  influenced  by  the  vapor  pressure.  The 
experimental  measurement  of  the  vapor  pressure  of  CL20 
should  be  performed.  Second,  the  conditions  that  appear 
to  be  most  viable  for  maximizing  the  solubility  of  CL20 
in  supercritical  carbon  dioxide  are  the  lowest  supercritical 
temperature  examined  (35°C)  and  the  highest  pressure. 

REFERENCES 

Ambrose,  D.A.  and  Walton,  J.,  1989,  Pure  &  Appl.  Chem.  61, 
1395-. 

Behrens,  R.,1990,  J.  Phys.  Chem.,  94(17),  6706-6718. 

Burov,  Y.M.,  Nazin,  G.M.,  Manelis,  G.B.,  1999,  Russ.  Chem. 

Bl.,  48(7),  1250-1254,  lav.  Akad.  NaukSer.  Khim  (1999). 
Dionne,  B.C.,  Rounbehler,  D.P.,  Achter,  E.K.,  Hobbs,  J.R.,  Fine, 
D.H.,  1986,  J.  Energetic  Materials,  4,  447-472. 

Edwards,  J.,  1953,  Trans  Faraday  Soc.,  49,  152. 

Joback,  K.G.,  Reid,  R.C.,  1987,  Chem.  Eng.  Comm.,  57,  233- 
243. 

Maksimov,  Y.Y.,  1992,  Russ.  J.  Phys.  Chem.,  66(2),  280-281. 
McHugh,  M.,  Paulaitis,  M.E.,  1980,  J.  Chem.  Eng.  Data,  25, 
326. 

Modell,  M.,  Reid,  R.C.,  1983,  Thermodynamics  and  Its 
Applications,  (2nd  ed„  Prentice  Hall  Publishing,  Englewood 
Cliffs,  NJ). 

Morris,  J.B.,  1998,  J.  Chem.  Eng.  Data,  43,  269-273. 

Nielsen,  A.T.,  Chaffin,  A.P.,  Christian,  S.L.,  Moore,  D.W., 
Nadler,  M.P.,  Nissan,  R.A.,  and  Vanderah,  D.J..,  1998, 
“Synthesis  of  Polyazapolycyclic  Caged  Polynitramines,” 
Tetrahedron  54,  11793-11812. 

Peng,  D.Y.  and  Robinson,  D.B.,  1976,  “  A  new  two  constant 
equation  of  state,”  Ind.  Eng.  Chem.  Fundam.  15,  59-64. 
Poling,  B.E.,  Prausnitz,  J.M.,  O’Connell,  J.P.,  2001,  The 
Properties  of  Gases  and  Liquids,  (5th  ed.,  McGraw-Hill, 
New  York,  NY). 

Sandler,  S.I.,  1989,  Chemical  and  Engineering 

Thermodynamics,  (2nd  ed.,  John  Wiley  Publishing  Co.,  New 
York,  NY). 

Schmitt,  W.J.,  1984,  The  Solubility  of  Monofunctional  Organic 
Compounds  in  Chemically  Diverse  Supercritical  Fluids, 
Ph.D.  Dissertation,  Massachusetts  Institute  of  Technology. 
Stein,  S.E.,  Brown,  R.L.,  1994,  J.  Chem.  Information  and 
Computer  Science,  34(3),  581-587. 

Suri,  S.C.,  Chapman,  R.D.,  1988,  Synthesis,  743-745. 
Tsekhanskaya,  Y.V.,  Iomtev,  M.B.,  Mushkina,  E.V.,  1964,  Russ. 
J.  Phys.  Chem.,  38(9),  1173. 

U.S.  Department  of  Labor,  Occupational  Safety  and  Health 
Administration,  “Cyclonite  (RDX)”, 

http://www.osha.gov/dts/sltc/methods/partial/pv2 1 35/pv2 1 3 
5.html>,  accessed  11/01/03. 


8 


